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Abstract 



Nuclear transparency is calculated for high-energy, semi-inclusive (e, e'p) 
reactions, by accounting for all orders of Glauber multiple-scattering and by 
using realistic finite-range pN interaction and (dynamically and statistically) 
correlated nuclear wave functions. The nuclear correlation effect is reduced 
due to the pN finite-range effect. The net effect is small, and depends sensi- 
tively on details of the nuclear correlations in finite nuclei, which are poorly 
known at present. 
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Typeset using REVT^ 



Perturbative quantum chromodynamics (pQCD) predicts the novel nuclear phenomenon 
of color transparency, in which a hadron produced by a hard scattering process could have an 
unusually long mean-free-path in nuclei |I| . Such a hadron would be small, color neutral, and 
suffer little interaction with other nucleons as it goes through the nucleus. The phenomenon 
has been observed in {p, 2p) and vector meson production ^ , while it is not seen in 
{e,e'p) pQCD calculations largely remain semi-phenomenological and have not yet 
yielded precise quantitative predictions. 

Low-energy hadrons can have long mean-free-paths in nuclei solely due to nuclear medium 
effects without color transparency It is imperative to establish nuclear medium effects 
reliably in these high-energy reactions, so as to identify the pQCD phenomenon explicitly. 
The medium effect that has been most noted is nuclear correlation, but calculations vary 
wildly from only a few percent increase to as much as a 20 to 30 percent increase in 
proton emission. 

The correlation effect results from an interplay between the nuclear reaction mechanism 
and nuclear structure. As such, both must be formulated and computed with equal rigor and 
care in order to unambiguously establish the significance of the effect. In this letter we report 
results of calculations of nuclear transparency in {e,e'p). We include all orders of Glauber 
multiple-scattering series with realistic finite-range pN interaction together with realistic 
nuclear correlations, and carry out numerical evaluation using Monte Carlo integration over 
all nucleon coordinates and spins. No previous work has been carried out with this rigor . 

The central quantity in this work is the nuclear transparency T, which corresponds to a 
ratio of the proton emission rate with and without the final-state interaction. T is given as 
the missing-momentum integration of the response function, 

= ~- < 0\O^ ■ ImG^ ■ O\0 > (1) 

TT 

for the ground-state nucleus |0 >. Here, ImG^ is the imaginary part of the Green's function 
due to proton emission, and 

ImG^ = (1 + GVi|^_i)JmGo(l + ^i.^-iG), (2) 
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where Gq and G are the Green's function for the proton in free space and in the nucleus, 
and Vi^A-i is the interaction between the knockout high-momentum proton (labeled as 1) 
and the rest of the A — 1 nucleons ||10|| . 

We apply the Glauber approximation for Green's functions, which consists of the eikonal 
approximation for the proton and the fixed-scatterer approximation for the target nucleons, 
and obtain 

/ ■ ■ ■ rfrAn/=2e-2^->^^('-^''-^)|^(ri, ■ ■ ■ , r^)!' 

I = 2 (3j 

Jdri ■ ■ ■rfr^|^f(ri, • • ■ ,r^)| 

with 

Xj(ri, r,) = -— / dz^V^t^iv^ - (4) 
Equation appeared in our previous work ||ll|] and and also in Ref. as an intermediate 



step to Glauber multiple-scattering series [12 



T can be cast in the well-known eikonal form expressed in terms of the one-body nuclear 
density, p(r): 

^ / c?rie3:p(-2/mx(ri))p(ri) 

Here, 21mxij\) is expressed in a cluster expansion as 

5:^^ J dr,(2i?er,-|r,f)p2(ri,r,)Mri) + ---, (6) 

where = 1 — exp(+2Xj(ri, rj)) is the Glauber A^A^ amplitude in the coordinate space. 

Nuclear correlation effects appear explicitly in the leading order through the two-body 
nuclear density p2, and spectator effects appear in the subleading order (not explicitly shown 
here) [0]. The commonly used eikonal expression is a consequence of invoking the zero-range 
approximation and setting Fj as a constant proportional to the A^A^ cross section, a: Either 

2Jmx(ri) ^ a / dz^g{Yx - r'i)p(r'i)b,^=bi (7) 

J 21 



or 



2/mx(ri) ^ a / dz[p{vW,^^^^ (8) 

with or without the correlation, repectively. Furthermore, many authors use the total cross 
section, atot-, as a. This corresponds to discarding the \T.j\^ term in Eq. and is not 



justified for a fully semi-inclusive process [|rT| |T3| since the term is of the leading order 
through unitarity |T^, |[Tl|. Note that the use of atot is also inappropriate in exclusive 
(e,e'p) reactions |n . 



VpN in Eq. (Q) is either the pp or pn potential, depending on whether the j-th nucleon is 
the proton or neutron. We determine VpN by applying the Abel transformation to invert the 
eikonal forms of the A^A^ amplitudes |T2[. The inversion is exact in the eikonal formalism 
for a local, spherical potential; we thus determine V^n better than the frequently used Born 
approximation. The A^A^ amplitudes are extracted from experimental data in Gaussian 



forms [1^. Because of the limited availability of pN scattering data, we also restrict the 
potentials to be spin-independent. The free-space pN interaction would be modified by 
nuclear medium corrections and by long formation lengths of excited nucleons. These are 
clearly of much interest, but we proceed by neglecting them in order to determine the 
correlation effects with minimal complication. 

We now describe our findings by taking the case of pi = 4.49 GeV/c for doubly closed 
shell nuclei, ^^O and ^°Ca. This choice is also made to minimize effects of other matters that 
might complicate the issue, such as open nuclear shells and strong momentum dependence 
of the pN amplitudes. Note that the pi value corresponds to the highest in the NE18 
(e, e'p) experiment . 

As Eq. (I) shows, exp(— 2/mx(r, 0)) contributes to T as a product with the correlated 
nuclear density. In Fig. 1, we illustrate exp{—2Imx{r,0)) calculated from Eq. using 
the Vpp inverted from the observed pp amplitude. Figure 1 shows that exp{—2Imx{r,0)) 
deviates from unity even for b and z > 1 fm. This feature is missing in the commonly used 
approximation of the zero-range pN interaction: The corresponding figure in the zero-range 
approximation would appear as a deep narrow valley following the z-axis starting from z = 0. 
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In order to evaluate T, we require the A-body nuclear density that is accurate in the 
length scale of several tenths of fm or smaller. Since such densities are not readily available 
for various nuclei, we construct them by following a simple scheme and compare with a 
variational Monte Carlo calculation of ^^O [|I^. We use a Jastrow-form of the wave function: 
a product of antisymmetrized Hartree-Fock (HF) single-particle wave function (for Skyrme 
force |T^) ^4$ and of pair wave functions for nuclear matter at various densities /'s, 



which describe the statistical and dynamical correlations, repectively. The A-body density 
|*(ri,---,r^)|2 is 

i<j 

±2(/e/. + SU^r + Gftftr - fr' " ^far' - 6ftr')] (9) 

for a pn pair (— ) and a pp or nn pair (+). Invoking the local density approximation, we 
evaluate the /'s as functions of |rj — rj|, for the (one-body) nuclear density that is generated 
from \A^\1 (spin-summed) at |rj + rj|/2. 

Here, we include the dynamical correlations only between nucleon pairs, neglecting corre- 
lation operator symmetrization and off-diagonal contributions. The terms depending linearly 
on the spin-dependent operators are discarded, and those depending quadratically are av- 
eraged. This simple scheme should be reasonable for closed-shell nuclei, especially because 
all noncentral /'s are much smaller than the central fc- The tensor correlations become 
stronger relative to fc at short distances (< 0.5 fm), affecting significantly the momentum 
distribution above the Fermi momentum ||20|, but this aspect is of little significance in our 



calculation using Eq. (^. 

With these quantities, T is computed by Monte Carlo integration of Eq. (|^). The integra- 
tions also sample the weight, i. e., the nuclear density. Figure 2 illustrates the two-particle 
distribution functions, pp„(|r|) and Ppp{\r\) for ^^O. For comparison, we also show the dis- 
tribution functions obtained using the without the /'s. Our pp„ and ppp have the same 
shapes as those by the variational Monte Carlo calculation with the peak being at the 



same location. But our peak heights are lower than theirs by about 40 %, and our ppn 



and ppp calculated from ^4$ alone are also lower than their mean-field results by a similar 
amount. Our one-body density changes little by the application of /'s and remains close to 
observation, while the one-body density by the variational calculation changes substantially. 

Figures 1 and 2 illustrate our main finding in this letter that the finite interaction-range 
effect reduces the nuclear correlation effect |]TU[. We elaborate on this in the followings. 

The values of T are tabulated in Table I, computed by the use of Eqs. (H) and (H) under 
various assumptions. T using Eqs. (0) and (H) is referred to as eikonal and is calculated 
with the one-body density obtained from Eq. (P). We show a)the most naive but popular 
calculation using atot and b) also with g{r) = 1 — exp(— (r/a)^) (a = 0.7 fm). T is increased 
from a) to b) by about 30 % and becomes accidentally rather close to that of our full result 
in 6), but this hides the essential physics that we are addressing. 

In 1) we show the eikonal result by the use of (Tinei- This corresponds to the leading- 
order [Eq. (H)] calculation of our full expression [Eq. (j^)] under the zero-range approximation 
without correlation, as it has been noted previously ||n| [0] |[13|. 1) is the one that should be 
compared with Monte Carlo results of the full expression, 2) - 6) labeled as GMS, Glauber 
multiple- scattering. 2) - 5) examine separately effects of the dynamical and statistical corre- 
lations. Once the finite-range interaction effect is included, we find a little variation among 
T's of 2) - 6). 6) is our full calculation. By comparing 1) and 6), we see that the net effect 
due to both nuclear correlation and finite-range interaction is of several %, above the eikonal 
T with ainei- This establishes our main finding. Note that the net effect becomes less in the 
heavier nucleus. 

In the rest of this Letter, we demonstrate that our main finding is affected little by the 
choice of the dynamical A^A^ correlations. As noted above, the peak heights of our pp„ and ppp 
differ from those by the variational calculation. While ours remain close to the observation 
within about 1 %, the RMS radius of the variational ^^O density is smaller than the observed 
by about 8 %. When we artificially enhanced our /'s so as to obtain the variational pp„ and 
Ppp, we observed our RMS radius decreased as much as the variational density did (though 
ours had a shallower central depression.) T computed with the enhanced /'s is 0.548 ± 
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0.003, which is smaller than our full calculation of 6), 0.577 ± 0.002. That is, the enhanced 
/'s causes nuclei to be less transparent. 

The reason for the effect is as follows: At a low nuclear density, fc around 1 - 2 fm becomes 
larger than the asymptotic value of unity, caused by the attractive, intermediate-range A^A^ 
interaction |T^. The large fc lifts the peaks in the two-particle distributions above the mean 
field results and increases the one-body density. Consequently, the knockout proton suffers 
more final-state interaction. The variational one-body density differs from observation, and 
our wave function is not constructed for the energy minimization. We conclude that the 
two-particle distributions are not reliably known, causing to be T uncertain, perhaps at 
about the 5 % level. 

In order to demonstarte this point, we also use a schematic fc, 1 — exp{—{r/af') (a = 
0.7 fm) by setting all other /'s zero. For this, Pp„ and Ppp are found to be always below 
the HF distributions, the RMS of the density increases from the HF value of 2.65 to 2.71 
fm, and T increases to 0.605 ± 0.003. What the correlation does in this case is merely to 
enlarge the nuclei and to make them more transparent. Note that the correlation effect 
usually discussed in the literature includes only this short-range repulsive aspect. 

In this work, we have examined the reaction that is fully semi-inclusive. That is, the 
Era integration is over the full range, and this permits the use of closure to simplify our 
calculation. A Monte Carlo simulation is underway to examine how realistic this procedure 
is by comparing with the NE18 experiment setup [^. Our calculations reported here are 



also for simplified kinematics without missing momentum (pm) cut. Our formalism allows 
Pm to be finite, though Eq. (|^) becomes more complicated, depending on one-particle off- 
diagonal nuclear density and the real part of VpN- Furthermore, our formalism is also 
applicable to inclusive (e, e'), involving the correlations differently from the semi-inclusive 
(e, e'p) case considered here. We will describe a detailed account of the present work and 
these considerations in forthcoming publications. 

In conclusion, we find that 1) the finite-range pN interaction reduces the effect of nuclear 
correlation, 2) the net effect is small and becomes smaller in heavier nuclei, and 3) the 



precise value of the transparency above the eikonal value depends sensitively on details of 
the dynamical correlations in finite nuclei, which are presently known rather poorly. 
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TABLES 

TABLE L Nuclear transparency of closed-shell nuclei at the knockout proton momentum 
4.49 GeV/c in various approximations and models: With or without pN finite-range interaction 
(FRI), dynamical correlation (DC), and statistical correlation (SC) in the eikonal and Glauber 
multiple-scattering (GMS) formalism. The numbers in parentheses denote equations used, g 
denotes the use of the simple form of the pair-distribution function (shown in text) and POD of a 
product of one-body densities generated by integarting the ^-body density indicated. 5) excludes 
DC among spectators (not involving the knockout proton.) Single-digit values in parentheses are 
statistical uncertainties (in the last digits) in Monte Carlo integrations. 





FRI 


DC 


SC 


Comments 


160 


40 Ca 


a) 








Eikonal [((0);CTtot] 


0.444 


0.335 


b) 








Eikonal[(||);(Ttot,fl'] 


0.529 


0.405 


1) 


Zero 


No 


No 


Eikonal[(0);(jj„ei] 


0.560 


0.446 


2) 


Finite 


No 


No 


GMS[POD;(|^^'P)] 


0.570(2) 


0.439(2) 


3) 


Finite 


Yes 


No 


GMS[POD;(|D] 


0.568(2) 


0.442(3) 


4) 


Finite 


No 


Yes 


GMS[[^^|2] 


0.579(3) 


0.447(3) 


5) 


Finite 


Yes 


Yes 


GMS[(|);Nospect] 


0.582(2) 


0.450(3) 


6) 


Finite 


Yes 


Yes 


GMS [(|); Full] 


0.577(2) 


0.447(3) 
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FIGURES 

FIG. 1. ej;p(— 2/mx(r, 0)) using the spin- independent, local Vpp extracted from experimentally 
determined pp scattering amplitude. The quantity illustrates attenuation of the knockout proton 
as it impinge on a proton located at the origin. The proton comes from the negative z, parallel to 
the z-axis. 

FIG. 2. The two-particle distributions of ^^O, ppn and ppp, calculated using the wave function 
of Eq. with and without pair correlation functions, shown by full curves and open circles, 
respectively. 
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